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Abstract 

We systematically develop the procedure of holographic renormal¬ 
ization for RG flows dual to asymptotically AdS domain walls. All 
divergences of the on-shell bulk action can be cancelled by adding co¬ 
variant local boundary counterterms determined by the near-boundary 
behavior of bulk fields. This procedure defines a renormalized action 
from which correlation functions are obtained by functional differenti¬ 
ation. The correlators are finite and well behaved at coincident points. 

Ward identities, corrected for anomalies, are satisfied. The correlators 
depend on parts of the solution of the bulk field equations which are 
not determined by near-boundary analysis. In principle a full non¬ 
linear solution is required, but one can solve linearized fluctuation 
equations to define a bulk-to-boundary propagator from which 2-point 
correlation functions are easily obtained. We carry out the procedure 
explicitly for two known RG flows obtained from the maximal gauged 
D=5 supergravity theory, obtaining new results on correlators of vec¬ 
tor currents and related scalar operators and giving further details on 
a recent analysis of the stress tensor sector. 
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1 Introduction 


The AdS/CFT correspondence S 1,1 has provided a way to obtain strong 
coupling results on conformal field theories from string theory and super¬ 
gravity. The techniques can be extended to describe non-conformal theories 
which are deformations of CFT's either by addition of relevant operators to 
the Lagrangian or by turning on VEV's for these operators. The gravity 
duals of such theories are domain wall solutions of (d + l)-dimensional bulk 
theories whose isometry group is the d -dimensional Poincare group. Many 
such classical solutions are known ! 1.1. 0. i, i. 0 |lTJ . Given such a 
solution one may address the problem of calculating correlation functions of 
operators in the deformed CFT. Much information about 2-point correla¬ 
tion functions has been obtained using a method which was formulated quite 
early in the game |3], |12[. But, as we discuss below, this method does not 
deal adequately with divergences and it is not straightforward to obtain cor¬ 
relators of the energy-momentum tensor and vector currents. The purpose 
of this paper is to apply the method of holographic renormalization |L3], [14 
to this problem. A finite renormalized action S ren can be constructed sys¬ 
tematically using this method. S ren is a functional of boundary data for the 
bulk fields which are the sources of operators in the field theory. Renor¬ 
malized correlation functions are then obtained by functional differentiation, 
and these correlators satisfy appropriate Ward identities including anomalies. 
The calculation of 1- point functions is particularly clear in this method. 

The domain wall spacetimes asymptotically approach AdSd+i geometries 
and the divergences which occur in the holographic computation of correla¬ 
tion functions are due to the infinite volume of these geometries. These IR 
divergences correspond to the UV divergences which occur in held theory 
when composite operators approach coincident points. 

We now outline how the IR divergences are resolved in the method of 
holographic renormalization. The first step is to express bulk fields as series 
expansions in the radial coordinate which is transverse to the boundary. This 
allows one to determine an asymptotic solution of the held equations given 
arbitrary Dirichlet boundary conditions. The solution is obtained by substi¬ 
tuting the series expansions in the nonlinear bulk held equations and solving 
term-by-term. This process is called near-boundary analysis. Roughly speak¬ 
ing it determines the hrst A — d/2 terms in the expansion where A is the 
dimension of the operator dual to a given bulk held. Beyond this point ad- 
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ditional information is required to determine a unique solution, as expected 
since we have second order field equations and only the Dirichlet condition 
at the AdS boundary has been specified so far. To specify a unique solu¬ 
tion we require that fluctuations about the domain wall background vanish 
in the deep interior. The coefficients determined by near-boundary analysis 
are local functions of the boundary data, while higher order terms can be 
non-local. 

The next step of the method is to construct S ren by a process of regulariza¬ 
tion and renormalization. Both the concepts and details of this construction 
are important. The bulk theory is regulated by cutting off the spacetime 
at a large finite value of the radial coordinate. The series solution is then 
inserted in the regulated classical action. One then observes that the on- 
shell action contains a finite number of terms which diverge as the cutoff is 
removed. These divergences involve only coefficients in the solution which 
are fixed by near-boundary analysis. They can be cancelled by adding coun¬ 
terterms to the action which are expressed as invariant local functionals of 
the induced metric and other fields at the boundary of the cutoff space-time. 
These fields depend locally on the Dirichlet data at the true boundary which 
is approached as the cutoff is removed. The entire process is similar to what 
is done in field theory, and there is an ambiguity of finite local counterterms 
which is just the usual scheme dependence in field theory. The sum of the 
regulated action plus counterterms is finite as the cutoff is removed, and S ren 
is defined by this limit. By virtue of its construction S ren is invariant under 
5D diffeomorphisms, except for the ones which generates Weyl transforma¬ 
tions of the boundary metric. The violation of Weyl invariance is manifested 
in logarithmically divergent counterterms, and conformal anomalies can be 
read directly from these terms. 

Finite correlation functions can then be obtained by functional differen¬ 
tiation of ^en with respect to sources. This computation, which involves a 
limiting procedure described in Sec 4 below, can be carried out in all gen¬ 
erality, and the correlation functions then involve the lowest order series 
coefficients which are not determined by near-boundary analysis. This is 
to be expected since correlation functions are non-local. By contrast, field 
theory UV divergences and anomalies are local, and as such one can obtain 
them by the near-boundary analysis. This also implies that Ward identities 
can be established by the near-boundary analysis. 

The domain-wall spacetimes that describe the RG-flows generically con- 
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tain naked singularities in the interior. Near the singularities the curvatures 
are large and the supergravity approximation breaks down. A full string the¬ 
ory treatment is presumably necessary to resolve the singularity. However, 
at least for quantities that do not depend crucially on IR properties, one 
should still be able to use supergravity. Taking this point of view, a very 
important issue is what boundary conditions are imposed at the singularity. 
In principle this information should come from string theory. Here we take 
as a working assumption that the fluctuations vanish at the singularity. 

Ideally one would like to have an exact solution of the nonlinear held 
equations with arbitrary Dirichlct boundary conditions and suitable behavior 
in the interior. The coefficients in the boundary expansion not determined 
by local analysis could then be extracted and n-point correlation functions 
computed. This problem is now intractable. However, it is not difficult 
to solve the fluctuation equations linearized about the background. The 
solutions are the extension to i?G-hows of bulk-to-boundary propagators in 
the pure AdS background f2|, and 2-point correlators are easily obtained 
from their boundary asymptotics. Higher point functions can be obtained 
by means of Witten diagrams, but we shall not pursue this here. 

T2| for calculating 2-point functions, one 


In the conventional method 


also imposes a cutoff at large radius and solves the linearized fluctuation 
equation with Dirichlet boundary condition at the cutoff (as opposed to the 
scaled Dirichlet problem on the full spacetime in holographic renormaliza¬ 
tion). The second variation of the cutoff on-shell action is then computed 
in momentum space yielding an expression containing singular powers of 
the cutoff times integer powers of p plus nonsingular terms in the cutoff 
which are non-analytic in p. The 2-point function is defined as the leading 
non-analytic term. The polynomial terms in p which are dropped are just 
contact 5-function terms in the position space correlator, which are scheme 
dependent in held theory and largely unphysical. The non-analytic term has 
an absorptive part in p which correctly gives the 2-point function for sep¬ 
arated points in x-space. Although this method is quite efficient, it is not 
fully satisfactory since it is not correct in general simply to drop divergent 
terms in correlation functions. Subtractions in different correlators should be 
mutually consistent, and this is guaranteed if the subtractions are made by 
local covariant counterterms. Further, to our knowledge, there is no system¬ 
atic method other than holographic renormalization to obtain correct finite 
values for 1-point functions. We hope to convince readers that holographic 
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renormalization is a formalism which resolves divergences in all correlators 
and makes Ward identities and their anomalies manifest ab initio. Further, 
the procedure need only by carried out once for each given bulk action with 
results applying to all classical solutions. Considerable insight into the rela¬ 
tion of holography and held theory is also exhibited, including the fact that 
the cancellation of UV divergences does not depend on the IR physics. 

In previous work [151 we have applied holographic renormalization to the 
computation of 2-point functions of the stress tensor and the operator dual 
to the scalar held in the domain wall background for two specihc i?G-hows, 
the GPPZ how [[7| describing Af = 4 SYM theory perturbed by an Af = 1 su¬ 
persymmetric mass term, and a Coulomb branch how [j|, |6j] describing a disc 
distribution of H3-branes. Essentially equivalent results were simultaneously 
presented in |T6[ with earlier progress in |17|. In this paper we present details 


of the near boundary analysis of the held equations which were omitted in 
|T5|j , and we extend the method to the computation of correlators of some 
vector currents. The emphasis here is on longitudinal components of the cur¬ 
rents, since the transverse 2-point functions can be obtained p~8| , pdf] by more 
conventional methods. However, we do summarize recent developments on 
the issue of Goldstone poles in transverse correlators. 

In Sections 2 and 3 we review the background solutions for the GPPZ and 
CB hows and the Lagrangians and equations of motion required to obtain 
the correlation functions we study. In Section 4 we describe the procedure of 
near-boundary analysis and how it leads to 1-point functions (with sources) 
that satisfy correct anomalous Ward identities and from which higher point 
correlators can be obtained. In Section 5 details of the near boundary analysis 
of the stress tensor sector of both hows are given. These results were already 
used in |L5| to obtain correlation functions. Sections 6,7,8 are devoted to the 
current sector. Near boundary analysis for vector helds, linear huctuations, 
and correlation functions are all discussed. In the Appendix a group theory 
lemma on the preserved havor symmetry on the Coulomb branch of Af = 4 
SYM theory is stated and proved. 


2 Holographic RG-flows 

In this section we review the method generally used to obtain supersymmetric 
domain wall background solutions of 5-dimensional supergravity theories. 
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The first step is to apply symmetry arguments |T0| to obtain a consistent 
truncation of the full bulk theory, usually D = 5, H = 8 gauged supergravity, 
to a small number of scalar fields interacting with gravity. In the cases we 
discuss a single scalar is sufficient and the truncated supergravity action is[] 


N 2 r 

S 0 = 7TT / 

2tt z Jm 


is + + v($) 



( 2 . 1 ) 


where /C is the trace of the second fundamental form. We work in Euclidean 
signature. We assume that the potential E(<f>) has a stationary point at 
<I> = 0. By a constant Weyl transformation m one may set E($=0) = —3, 
so that the action admits a pure AdS solution with scale L = 1 in string 
units. As discussed in JTH] , the Weyl rescaling produces the overall factor 


of N 2 /2it 2 in (|2.1|) , and the holographically computed correlation functions 
agree with those of the undeformed J\f = 4 SYM at short distance. 

The second step is to look for solutions of the Euler-Lagrange equations 
with d = 4 Poincare symmetry. The most general form is 


ds 2 = e 2A ^ Sijdx 1 dx^ + dr 2 , <f> = $(r) (2.2) 


Equivalent forms involving a different radial variable will be used later. With 
radial variable r the asymptotic boundary AdS region is r —> oo, and the 
scale factor is exponential in this region, i.e. exp(2A(r)) —> exp(2r). The 
asymptotic behavior of the scalar field partly determined by the mass or 
equivalently the dimension of the operator dual to it, the two being related 
by A = 2 +a/ 4 + m 2 , distinguishes solutions which are dual to flows governed 
by operator deformations from those triggered by VEV’s. 

Even within the ansatz (|2.2| ) it is generally difficult to solve the second 
order field equations of (|2.1| ). The following simplifying procedure is thus 
widely used. Namely if the potential V (<f>) is derivable from a superpotential 
W($), 


V 


1 dW 

2 ~d¥ 



(2.3) 


then any solution of the first order equations 


dA _ _2 W d&__dW_ 

_ dr _3 ’ dr <9<f> 

4 Our curvature conventions are as follows = d^VuK, X + r /ip A r j/K p — /j, 

Rfiv = R[i\v X - They differ by an overall sign from the conventions in |^o|, Q. 


(2.4) 

<-> v and 
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provides a domain wall solution for the action ( j2.1| ) [22| . p3| . These equations 
have been obtained from several standpoints, namely gravitational stability 
21, , fermion transformation rules in the (truncated) supergravity theory 

0], and the Hamilton-Jacobi framework [^I|. We will be applying (|2.3| ), 
in a supersymmetric context associated with Killing spinors, and the 


superpotential W will have a critical point at <3? = 0. 

We now review the application of this formalism to the two domain wall 
spacetimes for which we exemplify the formalism of holographic renormaliza¬ 
tion. It should be emphasized that the first order equations apply only to the 
domain wall background. The full second order equations of motion must be 
studied, via holographic renormalization, to obtain correlation functions of 
operators dual to fluctuations about the background. We begin this analysis 
in Section 4. 


2.1 Coulomb branch flow 

The first case we consider is the CB flow with n — 2, in the notation of |J. 
The behavior of the scalar field near the boundary (see below) shows that 
this solution corresponds to a VEV of a scalar operator of dimension 2. The 
VEV breaks R-symmetry from (6) to 50(2) x 50(4). The scalar that gets 
a VEV is the neutral singlet component of Tr (AWY- 7 ) in the decomposition 
20' (1, l)o © (1, l) +2 © (1, l) — 2 © (2, 2) +1 © (2, 2)_! © (3, 3)o under the 

50(4) x 50(2) subgroup of 50(6). Similarly, the vector currents decompose 
according to 15 —> (1,1) 0 © (1, 3) 0 © (3, 1) 0 © (2, 2) +1 © (2, 2)_i. Since the 
system describes a distribution of N D3 branes on a disk, the 50(2) symme¬ 
try is broken for any finite N. The corresponding held theory statement is 
that the Coulomb branch of J\f = 4 SYM with any finite dimensional gauge 
group G should not contain any configurations where the flavor group 50(6) 
is broken to a non-trivial direct product 50(n) x 50(6 — n ) with 2 < n < 5. 
We show in the appendix that this is indeed the case. 

As mentioned above, the how describes the theory at a point on the 
Coulomb branch of J\f = 4 SYM. The superpotential is given by 

_ 2 ® 1 44 > 

W(' h) --a v/e _ -ev'e. (2.5) 

Using (|2.3|), one can easily compute the potential. Near <3> = 0 it has the 
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expansion 


V(<f>) = -3 - 2<f> 2 + —^<f> 3 + C($ 4 ) (2.6) 

3y6 

which exhibits a tachyonic mass m 2 = —4, so that <f> is dual to an operator 
of scale dimension A = 2. 

The domain-wall solution is expressed implicitly in terms of a variable v. 
The solution and the relation between v and r are 

v = e VE *, e 2A = f— , d ^ = 2v 2 /\l-v) (2.7) 

1 — v dr 


The boundary is at v — l and the solution has a curvature singularity at 
v = 0. The parameter t is the radius of the disk of branes j6j in the 
10-dimensional “lift” of this solution. 

The scalar field vanishes at the boundary at the rate 


$ ~ 



1 

7e 


exp(—2r) 


( 2 . 8 ) 


For an operator of dimension A = 2, the dual bulk scalar approaches the 
boundary either at the rate r exp(—2r) associated with an operator deforma¬ 
tion of the Lagrangian or at the rate exp(—2r) associated with a VEV. We 
see that the VEV rate applies in this flow. 


2.2 GPPZ flow 

The second example we will deal with is the GPPZ flow [j7|. This solution 
corresponds to adding an operator of dimension 3 to the Lagrangian that 
gives a common mass to the three J\f = 1 chiral multiplets appearing in the 
decomposition of the J\f = 4 vector multiplet. The solution was proposed 
as the dual of pure Af = 1 SYM theory. Although it does not capture all 
of the expected features of the field theory, we will use it to illustrate the 
application of holographic renormalization to an operator deformation of a 
CFT A . 

The superpotential is given by 


W(4>) 


3 

4 


1 + cosh 



(2.9) 






Near $ = 0, the potential has an expansion, 

R($) = -3 - -<f> 2 - -<f> 4 + C>(<l> 6 ) (2.10) 

2 3 

The mass is m 2 = —3 indicating that the dual scalar has dimension A = 3. 
The domain-wall solution is given by 


$ 



1 + a/1 — u 
1 - a/1 — u 


e 2A = 


u 


1 -u 


( 2 . 11 ) 


where u — 1 — exp(—2r). The boundary is at u — 1 and the solution is 
singular at u = 0. Since <F ~ A/3exp(—r) near the boundary, we are dealing 
with an operator deformation by a dimension 3 operator, namely the top 
component of the superpotential AW = J2'j=i < h/- 


3 Lagrangian and Field Equations 


We will be interested in computing correlation functions in theories holo¬ 
graphically described by domain-wall backgrounds. This means that we need 
to solve the coupled system of equations that involves the corresponding bulk 
fields and then evaluate the on-shell action. For correlation functions that 
involve the stress tensor, certain R-symmetry currents and scalar operators, 
the bulk theory is that of gravity coupled to gauge and scalar fields. Since 
the present treatment is limited to 2-point functions, it is sufficient to keep 
only bilinear terms in the gauge fields. In this approximation the action for 
the SO(6) gauge fields reduces to a sum of uncoupled abelian sectors and 
it is convenient to use the Stiickelberg formalism involving gauge invariant 
combinations B M = A^ + d^a, with a bulk Goldstone fields. The correspond¬ 
ing held strengths are simply = d^B u — d v B^. (See |TR IHJ where the 


derivation of the vector sector Lagrangians from the full D = 5 gauged su¬ 
pergravity theory is discussed in detail.) The bulk Lagrangian that describes 
this system of fields is given by 


S = S 0 + / d 5 xVG 


IM 


-K($)F kX F, u G^G 


\v 


+ -M 2 ($)G^B^B u 


(3.1) 


where Sq is given in ( |2 . 1|) and A'($) and M 2 ($) are positive semi-definite 
functions of the active scalar $. We will only need their value along the solu¬ 
tions of the held equations. For the GPPZ how, up to a numerical rescaling, 


9 











the kinetic term is canonical and 


, r2 . . 2 / 2$\ (1 - u) 

M=smh Ur 

(3.2) 

For the CB flow 


Kso( 4) = ev / e # = u 2 / 3 

(3.3) 

K so{ 2) = e“75* = u- 4/3 

(3,4) 

K x = e~^ = u” 1/3 

(3.5) 

where X denotes the coset SO{&)/SO{4) x 50(2). 

The AdS masses vanish 


for the “unbroken” part of the R-symmetry group Mj 0 , 4 ) = Mj 0 , 2 ) = 0, 
while 


M\ = 4sinlr 


v / 6*h\ (1 — v)' 


The field equations that follow from this action are 


(3.6) 


^ ^ dV ldK 

□g® = ^ r KX r^ y 




<9$ 

-2 


4 <94> 

T — \ f ' TV T 

flU 2 '- T flV J- 1 ± 


G^ G ^F kX F uu + 


ldM 2 

- - G^BnBy , 

2 3$ ' ’ 


(3.7) 

(3.8) 

(3.9) 


where is the matter stress energy tensor and Tr T = G ,J ' v T jlu . 

As in past work on holographic renormalization it is most convenient to 
work in the coordinate system where the bulk metric takes the form 


ds 2 = G llu dx* J, dx l/ 


dp 2 
4 p 2 


H— gij(x, p)dx'xhd 


(3.10) 


where[] 


9ij { x ) P ) — 9(0)ij + 9{2)ijP + p 2 \g(4)ij + ^1(4 )ij log p + h 2 (4)ij{\og p) 2 ) + ... (3.11) 

5 In the Fefferman-Graham framework the most general expansion may contain half¬ 
integrals powers of p, or integral powers of the coordinate U, p = U 2 (the p variable was 
introduced in |H|). In the case of pure gravity and all the cases we consider here, the 
coefficients with odd powers of U can be shown to vanish (except the g(d)ij coefficient in 
d + 1 dimensions with d odd). See however p(i| for an example in d + 1 = 3 that involves 
a near-boundary expansion with half integral powers of p. 
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Any asymptotically AdS metric can be brought to be of this form near the 
boundary [27]. The radial variables p and r are related by p — exp(—2r). 
The boundary is located at p = 0 and the regularized action will be defined 
by restricting to p > e. 

In the coordinate system ( |3.1C|) , the scalar held equation reads 


1 OT / 

4P 2 4" - 4p4'(l - -p(logg)') + pD.t ~ — = 
P 2 dK, m ntic v , o_ kl c „ , , p3M : 


(3.12) 


4 a4> [g m "g u F mk F nl + 8 

where prime indicates derivative with respect to p, g = det gij(x, p), and 
is the scalar Laplacian in the metric g t -j(x, p). 

Keeping terms which are at most quadratic in the gauge fields and drop¬ 
ping R-symmetry indices, the vector held equations (|3.7| ) read 


di{\fgKg lJ F ]p ) = M 2 ^-B p 

Jjr, kl 


diiy/gKg^Fji) + 4pd p (^Kg kl F pl ) = M 2 —g kl B, 




ki 3 


p 


Taking the divergence of fl3.14|) we obtain 

9i(V9Af 2 9 «B,-) = 0 


(3.13) 

(3.14) 

(3.15) 


The Einstein equations in the coordinate system (|3.10|) are given by 
P - 2(sW)ii + Tr 0TV ) g'ij } + Rij[g} + 2 g'ij - Tr (< g~ l g') g tj = 


-2 + — [K(<h) - V(0)] gij + M BiBj + pK[g kl F lk F 3l + 4 P F ip F jp \ 


--pK gij [g mn g kl F mk F nl + 8 pg kl F kp F lp ] 


V,;Tr (g- l g') - = -2 (d&d p $ + M 2 B t B p + P Kg kl F lk F pl ) 


(3.16) 

(3.17) 


Tr Or 1 /) - -Tr (<ry<rV) = -2 <9 p 4>«9 p 4> + —JK(4>) - K(0)] + M B p B 


6 p 2 


%,Ktl u F t „F,„ - ^I<g mn g a F mk F nl 


(3.18) 


We are now ready to study the near-boundary solution of the above equa¬ 
tions. 
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4 Near-boundary analysis: generalities 


In this section we present a general discussion of the near-boundary anal¬ 
ysis. The emphasis is on the general structure rather than the details. In 
subsequent sections we specialize the analysis to the two specific flows we 
described in section 2. Although it was useful to distinguish between back¬ 
ground fields and fluctuations in the previous discussion, the near boundary 
analysis applies to the full bulk fields of the action fl3.1|) and must be later 
linearized with respect to the background. 


4.1 Asymptotic solutions 

The first step in the analysis is to find asymptotic solutions of the bulk held 
equations. To simplify the notation we suppress all spacetime and internal 
indices and collectively denote bulk fields by J-(x, p). Near the boundary 
each held has an asymptotic expansion of the form 

F(x } p) = p m (yf(p)(x) + f(2)(x)p + ••• + P n (f( 2 n)(x) + log p f( 2 n){x) + ...) + ...) 

(4.1) 

The held equations are second order differential equations in p, so there are 
two independent solutions. Their asymptotic behaviors are p m and p m+n , 
respectively. For bosonic helds in domain wall solutions dual to hows driven 
by protected operators deforming J\f = 4 SYM theory, n and 2m are non¬ 
negative integers. The boundary held /( 0 ) that multiplies the leading be¬ 
havior, p m , is interpreted as the source for the dual operator. In the near¬ 
boundary analysis one solves the held equations iteratively by treating the 
p-variable as a small parameter. This yields algebraic equations for fyk), 
k < n, that uniquely determine f^k) hi terms of f(o)(x) and derivatives up 
to order 2k. These equations leave f( 2 n)(x) undetermined. This was to be 
expected: the coefficient f( 2 n){x) is the Dirichlet boundary condition for a 
solution which is linearly independent from the one that starts as p m . As we 
will shortly see, f^n) is related to the VEV of the corresponding operator. 
The logarithmic term in (|4.1| ) is necessary in order to obtain a solution. It 


is related to conformal anomalies of the dual theory, and it is also hxed in 
terms of /( 0 )(x). 

As we have remarked, any solution which is asymptotically AdS can be 
brought to the p-coordinate system near the boundary. There are, however, 
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bulk diffeomorphisms that preserve the form of the coordinate system ( |3.10|) , 
but induce a conformal transformation at the boundary |28|. These trans¬ 
formations are related to the conformal anomaly [|T3J and thus with local 
RG equations on the field theory side, as discussed in a preliminary way in 
)30| , [T5H . As we will discuss, correlation functions are expressed in terms of 
certain coefficients in ( |4.1|) . It follows that the RG equations are encoded 
in the transformations of the coefficients under the bulk diffeomorphisms 
Such transformations for various bulk fields have been given in 


given m 


[28 


n n, m m . 


4.2 Counterterms 


The asymptotic solution obtained above can be inserted in the regulated 
action. A finite number of terms which diverge as e —» 0 can be isolated, so 
that the on-shell action takes the form 

*Sre g [/(o); e ] = [ d 4 x^/ c/( 0 )[e l 'a( 0 )+e < ' !/+1 ^a( 2 ) + ... — loge a^ + C^e 0 )] (4.2) 

J p=e 

where v is a positive number that only depends on the scale dimension of 
the dual operator and a^k) are local functions of the source(s) /(o). The 
counterterm action is defined as 


V [v.r .t = -divergent terms of 5 reg [/(o);e] 


(4.3) 


where divergent terms are expressed in terms of the fields jF(x, e) ‘living’ at 
the regulated surface p = e and the induced metric there, 7 ^ = gij(x,e)/e. 
This is required for covariance and entails an “inversion” of the expansions 
(|3.11|) ,(|4T|) up to the required order. 

To obtain the renormalized action we first define a subtracted action at 
the cutoff 

Ssubt-T 7 ^, e); e] = S ieg [f( 0 y, e] + Sa[F(x, e); e]. (4.4) 

The subtracted action has a finite limit as e —> 0, and the renormalized action 
is a functional of the sources defined by this limit, i.e. 


SrenL/io)] = hmS' sub [J r ; e] 

e—>0 


(4.5) 


The distinction between S su b and S Ten is needed because, as described in 
Sec |4.3|, the variations required to obtain correlation functions are performed 
before the limit e —> 0 is taken. 
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The procedure above amounts to a “minimal” scheme in which the diver¬ 
gences of S'reg are subtracted. As in standard quantum field theory, one still 
has the freedom to add finite invariant counterterms. These correspond to a 
change of scheme. For example, such finite counterterms may be needed in 


order to restore some symmetry ( e.g . supersymmetry [jT5 


Given a bulk action there is a universal set of counterterms that makes 
the on-shell action finite for any solution of the bulk field equations with the 
same Dirichlet boundary data. The counterterms are different for different 
bulk actions, for example for different potentials V($). 

The counterterm action can be decomposed into a set of counterterms 
associated to each field plus cross-terms. The former are universal in the 
sense that they are present everytime the bulk action contains a given bulk 
field, and the latter vanish if all bulk fields but one are set to their background 
value. This structure follows from the fact that the sources are independent. 
If one wants to study the theory with only one source turned one, one can 
either set all other sources to zero from the beginning and then compute 
counterterms or first compute counterterms with other sources non-vanishing 
and then set the other sources to zero. Quite obviously the two procedures 
should lead to the same result. In actual computations however, as the ones 
presented later on, this structure of the counterterm action comes about 
in a very non-trivial way. The dependence of the divergences in S reg on a 
given source changes when different number bulk fields are involved in the 
computation. However, the relation between the sources and the induced 
fields also changes when more bulk fields are involved in the computation. 
The combined effect is that the counterterms take the form discussed above. 

If the bulk held equations are coupled then one has to study the asymp¬ 
totic solutions of the coupled system of the equations in order to obtain the 
correct set of counterterms. In some cases, however, the near-boundary equa¬ 
tions decouple up to the order that determine coefficients that contribute to 
divergences. In these cases it is sufficient to study each field separately in a 
fixed background. This may simplify computations if one is only interested 
in computing a subset of all possible correlation functions. An example of 
this is the case scalars coupled to gravity in the Coulomb branch flow dis¬ 
cussed in section 5.1. If one only wants to compute correlation functions of 
the operator dual to the active scalar, the near-boundary analysis in section 
5.1.1 is sufficient. Only if the one wants to study correlation that involve the 
stress energy tensor as well one needs to analyze the full set of gravity-scalar 
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equations. This should be contrasted with the case of the scalars coupled to 
gravity in the GPPZ flow, discussed in section 5.2. In that case (which is the 
generic one) the gravity-scalar equations do not decouple in near-boundary 
approximation, so even when one is interested in only computing correlation 
functions of the operators dual to the active scalar, one still has to analyze 
the full set of the gravity-scalar equations. 


4.3 One-point functions 


Having obtained the renormalized on-shell action one can compute correla¬ 
tion functions by functionally differentiating S ren with respect to the sources. 
After the renormalization procedure discussed above variation of (|3.1|) reads 


fiS ien \g(p)ij, 0(o )j -4(o)*, °(o)] — / d 4 Xy/g(o)[-^(Tij)5g^ + (O$)50( O ) 


+ ( J l )5A(Q)i + (O a )5a(o)] 


(4.6) 


where g^, A( 0 )j, </>(o), U(o) are sources for the dual operators and appear as the 
leading coefficients in the near boundary expansions of the bulk metric G^, 
gauge field A^, active scalar and Stiickelberg field a, respectively. 

The expectation value of any scalar operator, such as the dual to is 
defined by 

<0 4 > = (4.7) 

VW) 00 ( 0 ) 

It can be computed by rewriting it in terms of the fields living at the regulated 
boundary^ 

' 1 1 SSsub \ 


( 0 $) = lim 


(4,8) 


e A / 2 yy 5$>(x, e)) 

where yy (x) = g VJ (x, e)/e is the induced metric on the boundary and 7 = 
det( 77 ). This expression allows us to compute the 1-point function in terms 
of (undetermined) coefficients in the asymptotic expansion of the bulk held 
$. 

The expectation value of the stress-energy tensor of the dual theory is 
given by 


/rr v 2 T5 ren 2 SS m b 

ij, = — - —^ = imi , „ ... -= hm 

VmSgl ;> -» 


Tab] 


(4.9) 


6 For scalars of dimension 2, an additional loge is needed in this formula, see (5.31). 
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where T tJ [ 7 ] is the stress-energy tensor of the theory at p = e. From the 
gravitational point of view this is the Brown-York stress energy tensor 
supplemented by appropriated counterterms contributions E3 


31 


Tub] = T’P + T§, (4,10) 

T-j S comes from the regulated bulk action and it is equal to (see (3.6) in B) 

Tir\l] = -1(K«-K 7 y) (4.11) 

1 3 

= ~2 (- d e9ij(x, e) + 9 ij(x, e) Tr [g~ l (x, e)d e g(x, e)] + -g^x, e)) 

where /Qj is the extrinsic curvature tensor. TF is the contribution due to 
the counterterms. 

Similarly, the 1-point function of R-symmetry currents reads 


(j) 


1 SS Ten _ ^ ( 1 1 5S snh \ 
^/W) 5A m ^ \e 2 v'T 5Ai(x , e) ) 


(4.12) 


The 1-point functions (C?$), (T^) and (J*) depend on all sources. Field 
theory vacuum expectation values can be obtained by setting the sources to 
zero. The same remarks hold for (O a ). However, it is somewhat formal to 
speak of a held theory operator O a dual to a Stuckelberg scalar held. It is 
better to convert to a linear formulation of the corresponding symmetry, and 
this is done in the examples discussed in Sec. |] 


4.4 Two-point functions 


In order to obtain higher point functions we need to functionally differentiate 
1 -point functions with respect to the sources. 

For example one can define and compute the connected 2-point functions 
of the stress tensor as 


i T ij(x)Tki(y)) = 


8(Tv(x)) 


fgwiy) ){y) ’ 

and the connected 2 -point function of the currents as 

1 S(Ji(x)) 


( Ji{x)J k (y )) = 


9(p){y) SA (o)(y) 


(4.13) 


(4.14) 
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In terms of the projectors 


PiPj rprpl 1 

0{j — i fcZ 7\\^ik^jl "F '^il'^jk) T^ij^kl (4-15) 

P Z o 

the (Fourier transform of the) former decomposes into 

(- Tij(p)T k i(-p )) = Uj^ kl A(p 2 ) + iTijiruBip 2 ) (4.16) 

while the latter into 

( Ji{p)Jk{~P )) = KikC(p 2 ) + ^ D(p 2 ) (4.17) 

pZ 

Supersymmetry Ward identities are expected to relate A to C and B to D , 
but we shall not derive these relations here. 


4.5 Ward identities 


Bulk gauge fields couple to boundary symmetry currents. It follows that 
bulk gauge invariance translates into Ward identities of the boundary quan¬ 
tum field theory. We derive here the Ward identities that the 1-point func¬ 
tions should satisfy. As we will explicitly show in the next few sections, the 
holographically computed 1-point functions do satisfy these Ward identities 
including anomalies. 

Using ( |4.6| ), invariance of (34) under diffeomorphisms, 


sg%) = -W + v j e ), <ty ( o) = rvi0 (o) , 

Sa ( o) = CVjO(o), 5A( 0 )j = ^VjA^ + Vj£-M( 0 )j, (4.18) 

implies the Ward identity for the conservation of the stress tensor)]. 

= —(O$)Vj0( o) — (O a )\7 jCi( o) — F(opj(J l ) + A(o)jV;(J') (4.19) 


where F^j is the field strength of A( 0 )j. 

7 Notice that these Ward identities are valid in the presence of sources. In particular, 
(Tij) depends on sources. In field theory one usually expresses the Ward identities in terms 
of the stress energy tensor with the sources set equal to zero, (Tij)QFT = (Uj) |sources=o- 
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Invariance under Weyl transformations, 

SgV = -2ag l \ (50 ( o) = -(4 - A)cr0 (O ) 

<54L(o )i = — <jy4( 0 )j, ^fl(o) = —<70(0) 

leads to the conformal Ward identity 

{Tt) = -(4 - A)0 ( o)(O„) - a (0) (O a ) - (J l )A m + A. 

where we have allowed for the conformal anomaly A. As explained in A 
is obtained directly from the logarithmic counterterm of the the bulk action. 

In the case of the GPPZ flow, the bulk Stiickelberg gauge invariance 
implies a corresponding Ward identity. We will find convenient in section 6.2 
to use normalizations such that B M — 3(5^002. Then the bulk gauge 

invariance implies for the sources 


(4.20) 


(4.21) 


‘( 0 ) 


It follows 


— V’A, 1(1(0) — gA 

(4.22) 


(4.23) 


A more detailed derivation of this Ward identity is given in section 6.2. 

In the Coulomb branch case, the renormalized on-shell is not invariant un¬ 
der transformations analogous to (|4.22|) because of the presence of Goldstone 
poles. We refer to section 6.1 for a detailed discussion of this point. 


5 Gravitational sector 

5.1 Near-boundary analysis for the CB flow 


The coordinate transformation that brings the domain-wall metric ( |2.2| ) to 
the coordinate system (|3.1C|) admits an expansion as [] 

1 - v = f p - -A p 2 + C>(p 3 ) (5.24) 


8 The exact coordinate transform is 


l 2 p = — 303 exp 


-A- - V3 tan -1 -L(i + 2 W 3 ) 


L03 


03 


(' 


, 1/3 _ 


1 ) 


01 + H 1 / 3 + V 2 / 3 
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In these coordinates the solution for $ and A, is given by 

i = ±(-f p+ by + o(p 3 )), e 2J = 1(1 - T t‘p 2 + 0(p 3 )) (5.25) 

v6 0 p 18 

The active scalar has dimension A = d/2 = 2 and can be decomposed as 

<h(x, p) = p log p <j>{x, p) + p0(x, p) (5.26) 


5.1.1 Scalars in a fixed background 


Let us first study the case of a scalar field in a given background of the form 
(13.10 ). As we will see in the next section, the backreaction of the scalar 
field to the gravitational field equations does not affect the coefficient of the 
metric that determine the divergences of the on-shell action, and vice versa. 
This means that for the purpose of the near-boundary analysis the scalar 
equation decouples from Einstein’s equations, so we can consistently study 
scalars on a fixed background. This means that the results of this section 
are also part of the analysis of the next section. 

The relevant field equation to solve is (|3.8Q with the potential in 
We look for a solution of the form (|5.26|) where 


(j){x,p) = 0(0) +0(2)P + plogp0( 2 ) + ..., (5.27) 

00*5 p) — 0(0) + 0(2 )P + ••• 


0(o) is the source for a composite operator of dimension A = 2, and 0(o) is 
proportional to the VEV of this operator (as we will shortly derive). 

The field equation ( |3.8|) can be solved order by order in the p variable. 
One obtains 

1/ 2 \4 9 1~ 

0(2) = -- (0(0)0(0) + 20(O) J R[5 , (O)]J - yg(0(o) - 20(0)0(0)) 

0(2) = ^D(O)0(O) + g-R[5 l (O)](0(O) + 0(0)) + 8(0(2) + 0(2))j 

0(2) = ^0(o) (5.28) 

As anticipated, 0( O ) is not determined by these equations. 
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Following the procedure outlined above, the regularized action is given 


by 


reg 



/ 


d d Xy/g( 0 ) log 2 e0 2 O ) + loge(0 2 O ) + 20( O )0(o)) + O(e 0 ) 


/ 



(5.29) 


Notice that the counterterm proportional to 1/loge is not vanishing in the 
UV, i.e. in the e —> 0 limit. This is so because it is <p^(x) rather than 
<F(a;, e) that is kept fixed in this limit. By inserting the expansion of y(x, e) 
and <F(;r, e) into the last expression in (|5.29|) one easily verifies that the infinite 
terms reproduce that of the second expression in (|5.29|) . The latter, however, 
differs from the former in its finite part. 

The renormalized action is given by 



$ 2 (x, e) 
loge 


(5.30) 


It is important that the counterterms are expressed in terms of field living 
on the regulating hypersurface. Otherwise, the subtraction would not be 
covariant. Although we have presented the derivation for D = 4 + 1 bulk 
dimensions, the final result is valid for operators of dimension A = d/2 in 
(bulk) dimension D = d + 1. The same holds for most of the following results 
in this subsection. 

Let us now compute the 1-point function, 



(5.31) 


where the first equality is the definition of the 1-point function and in the 
second equality we express things in terms of the regulated theory at p — e. 
Using Q5.3QD we obtain 



(5.32) 


20 




















where the last two terms is the contribution due to counterterms. Inserting 
in ( |5.31| ) one obtains 


( O $) — 20(o) 


(5.33) 


This is the relation advocated in [[h|] 9 . From this relation we conclude that 
the solution in (|5.25|) indeed describes a different vacuum state rather than 
a deformation of the original theory. Setting 0( O ) = —1/\/6, its background 
value, one finds the physical (0$) = —2/\/Q of this Coulomb deformation 
of the J\f = 4 theory. We view the derivation of vevs as a clean and elegant 
application of holographic renormalization. Had one differentiated the reg¬ 
ulated action (5.33) wrt 0( O ) and then discarded infinities, one would have 
gotten an incorrect result. This is so because the covariant counterterms 
include finite parts, as already remarked. 


5.1.2 Scalars coupled to gravity 


We now turn to study the backreaction of the scalars on the gravitational 
sector. To solve the equations we insert the asymptotic expansions in Ein¬ 
stein’s equations and solve them order by order in the p-variable. The most 
efficient way of performing this computation is to first differentiate the equa¬ 
tions with respect to p and then set p equal to zero. For instance, if we want 
the equations that appear at order p 2 , we may differentiate the equations 
twice with respect to p and then set p = 0. The resulting equations are alge¬ 
braic and can be readily solved to obtain the coefficients in the asymptotic 
expansion of the bulk fields. 

To leading order, one sets p = 0 in Einstein’s equations, solves them and 
finds that the scalars do not contribute and gm obtains the same value as 
for the case of pure gravity, 


9(2)ij 



(5.34) 


Here and henceforth curvature tensors and index contractions pertain to 
9(0)ij- 

The next order equations determine hi( 4 ),h 2 ( 4 ), the trace of p( 4 ) and the 
covariant divergence of p( 4 ), 


^ 1(4 )ij — h (4 )ij — -0(o)0(o)S , (o)«p h 2 (4 )ij — —7^4>lo)9(p)ij (5.35) 

9 For A = 2 scalar operators, the factor 2A — 4 is replaced by 2. 
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Tr 9 ( 4 ) = ^Tr g\ 2) - |(^ 0) + 2^ 0) ), (5.36) 

^ j 9m = V J (“[Tr 9(2) - ( Tr 9(2)?} 9(o)ij + \{g 2 ( 2 ))ij ~ \ 9(2)ij Tr g (2 ) 
_ g( 2 0(o) + 0(o))0(0 )ifj - (0(o)Vj0( O ) - 0(0)Vj0( O )) (5.37) 

where 

^(4 )ij = g RikjlR + — VjVj-R — — — A-Rjj 

+ (—V 2 A + —7? 2 - —i2 fcI i2 w ) = -T“ (5.38) 

l 96 96 32 u ’ 2 v ; 

and T)“ is the stress energy tensor derived from the action 

Sa f d XyJQ(0)-A.grav ■ (5.39) 


Ag r av is the gravitational trace anomaly given in (|5.48|) . The solution of the 
scalar held equation is described in the previous subsection and the results 
given ( |5.261 )- (|5.271) - (|5.28|) carry over. 

Having obtained the asymptotic solution we proceed to calculate the reg¬ 
ularized action. 


Sreg = 


[ d 5 xVG \]r + + H($)l -f- [ (5.40) 

2 p>e [4 2 J 2 Jp=e 


where 7 is the induced metric at p — e. A direct computation yields 


S reg = 


\J9( x , P) - 


/ d P dx 93 

lp>e 2 p 6 

3 

2 e 2 


— V 
3 . 


2 Jp=e 
2 


d 4 Xy/^JC 


(5.41) 


= J ~7^ ~ lo geg[(Tr g (2) ? - Tr gf 2) ] - log 


In the hrst equality we used Einstein’s equations to eliminate the scalar 
curvature from the action, and in the second equality we have inserted the 
asymptotic solution we just obtained. The logarithmic divergence agrees 
exactly with the one determined in (|5.29|) , and the gravitational one with 
the results of JT3| (as they should). 
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In order to compute the renormalized action ( |4.4| ), the divergent terms 
(as e —> 0) in (|5.40|) are to be cancelled by adding counterterms S ct = — Sdiv- 
When expressed in terms of the fields living at p — e, the result is 


Sot = 


I p=e 


d “ x ^ (1 ~ + ® 2(l ’ e >+ t ' 2(X ’ Cl 


fog t / 

(S-42) 

To derive this result one needs to invert the relations between the induced 
fields 7 ij(x,e), and the sources (?(o)p- and 0( O ). Details on how to do 

this are given in appendix B of 114]. The simplest way to check ( |5.42 ) is 


to to insert the expansions of 7 and $ and show that the divergent terms 
agree with the ones in ( |5.41|) . Notice that the log 2 e divergence due to the 
$ 2 term cancels against a similar divergence originating from ^ 7 . The first 
three counterterms in ( |5.42|) are the gravitational counterterms (see (B.4) of 
PH). Notice that the matter dependent infinities in ( 5.41|) are different from 
the ones in (|5.29|) . Nevertheless, the required matter counterterms coincide 
with the ones we computed in the previous section. The difference between 
the infinite parts is exactly compensated by the backreaction of the scalars 
to the metric. 

We can now proceed to compute the renormalized stress tensor. The 
contribution due to counterterms is given by 

<& 2 (x, e)' 


^ 




7 ij $ (x, e) + 


loge 


Combining (|4.11|) and (|5.43|) we obtain 


2 lo ^ T Z 

(5.43) 


(Tij) = -- lirn 


2 


-(-9(2)ij + 9{0)ij^ 9 ( 2 ) + 2 R ij - -^9(p)ijR) 


log e 


-2/i 2 (4)p - - 4>\o)9{o)ij 


+ l°ge ( —2(/ii( 4 )jj + h 2 (4)ij ) + — -(20( O )0( O ) + 0( O ))fi , (o)p 

~^9(4)ij — hi(4)jj + ( 7 ( 2 )ijTr g^ 2 ) ~ xfi , (o)qTr g ( 2 ) 


(5.44) 


kl 


■^(RikR k j — ZRikjiR 

O 

-~^9{2)ijR + g9{0)ij(RklR 


- VjV 7 f? + W R 


v 


6 


^ 2 Rg{ 0 ) 


ijj 


kl 


-R 2 ) 
6 ’ 


O ( 2 0(O) + 0(0) - 20(o)0(o))s , (o) 
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where the curvature tensor is that of g( 0 ). Using (|5.34| )- (|5.35|) we hnd that 
the 1/e, loge and log 2 e divergences cancel. The finite part is the expectation 
value of the stress tensor in the dual QFT. It is equal to 


(Ti S ) 


am + g[Tt 9(2) - ( Tr 9(2)) 2 ] 9(0)0 - |(9(: 2 ))« (5-45) 

1 x _ 2 3 

4 g( 2 )ij Tr p( 2 ) + -(0(o) - 30(o)0(o) )d(o)ij + g0(o)5 , (o)ij + 


1 T“ the contribution in the boundary 


Notice that because h^j is equal to 2 ^ l3 
stress energy tensor proportional to h^j is scheme dependent. Adding a 
local finite counterterm proportional to the gravitational trace anomaly will 
change the coefficient of this term. Similar remarks apply to the penulti¬ 
mate term: it is proportional to the stress energy tensor derived from an 
action equal to the matter conformal anomaly (|5.49|) . One may remove such 
contributions from the boundary stress energy tensor by a choice of scheme. 

Even though ^( 4 )^ is not fully determined by the bulk held equations, 
its divergence and trace are, and this is sufficient in order to compute the 
trace and divergence of (T i3 ) given in (|5.45|) . A direct computation using the 
asymptotic solution yields, 


(V) 


-(O$)V i 0 ( 0 ) (5.46) 

-20( O )(O*) + ~ ^R 2 ) + 20(o) (5.47) 


i.e. (Tij) correctly satisfies the diffeomorphism and trace Ward identities. 
The last two terms in (|5.47|) are what we called A in (|5.36| ). The second 
term 


Agrav = — (RijR lj - -R 2 ) 


(5.48) 


is the holographic gravitational conformal anomaly |T3| and the last term 


A scal = 20 2 o) (5.49) 

is the conformal anomaly due to matter |34| . The coefficients in both of them 
are known not to renormalize, and indeed we obtain the correct value. The 
Ward identities and the anomalies are important checks of the intermediate 
computations and of the consistency of the formalism. 
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5.2 Near-boundary analysis for the GPPZ flow 

In the coordinate system (|3.10| ), in which p — 1 — u, the functions 4> and A 
that determine the GPPZ domain wall solution expand as 


4> = / ,2 (V3 + -U+0(p 2 )), e 2A = --l 
V3 P 

Similarly, the expansion of the superpotential reads 

w = -\-\p + o(p 2 ) 


(5.50) 


(5.51) 


In our units, the quadratic term in the potential has m 2 = —3 implying that 
the operator dual to $ has scaling dimension A = 3. 

The gauge kinetic function of the 17(1) graviphoton is K = 3 and its 
effective mass term |[I8H is 


M 2 = -2A" = 4 


P 


(1-P) 5 


4[p + 2p 2 + 0(p 3 )\. 


(5.52) 


5.2.1 Scalars coupled to gravity 

In this section we perform the asymptotic analysis of the scalar-graviton 
system. These results have been used in []T5| without a proof. 

The asymptotic expansion of the metric is given in ( |3.10| ). The expansion 
of the scalar held reads 


<h(x, p) = p 1/2 (V>( 0 ) (x) + p(/>( 2 ) (x) + P log P0( 2 ) (x) + ■ ■ •) (5.53) 


Plugging ( |5.53|) into the scalar held equation one gets 


^(2) = -^ V2 0(O) + g-R[0(O)]0(O) 


(5.54) 


while the dependence of 0( 2 ) on the sources 0( O ) and g ( 0 ) remains undeter¬ 
mined. 

Einstein’s equations to the lowest order yield 

9{2)ij = ^ (-fyo )ij ~ q R (o) 9{ o)ij) ~ >(o)9(o)ij ■ (5.55) 
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From the computation of divergences in pure gravity we know that gravity 
divergences depend on g^ 2 ). Since there is a backreaction of the scalar to gp), 
these divergent terms lead also to scalar dependent divergences. These terms 
would be missed in a fixed background computation where g ( 2 ) is fixed and 
independent of the scalar fields. This is an example of the general discussion 
at the end of section 4.2. In cases like this, one needs to solve the coupled 
system of equations. Notice that in the Coulomb branch case the scalar field 
did not backreact on any of the metric coefficient that contribute to diver¬ 
gences, so it was consistent to consider the problem on a fixed gravitational 
background. 

The non logarithmic terms in Einstein equations yield 


Tl'g(4) — —20(0)0(2) — -0(O)Do0(O) — —R(p 2 ( 0 ) 


+—{ R ij RlJ - - R 2 ) + Q0(O) 


16 9 9 

V J c/( 4 )p- = V J (--[Tr gf 2) - (Tr g {2) ) 2 } g (0)ij + ^(c/ ( 2 2) )p- 


4 

+9(o)ij 


9 ( 2 )ij Tr g( 2 ) 

1 


'^1(4 )ij o Vj0(o) Vj0(o) 


(V0(0)) - 0(0) (0(2) + 0(2)) 


(5.56) 


-[-2(0(2) + 0(2)) + g0(o)]Vi0( O ) 


(5.57) 


The other coefficient tensor hu^j in the expansion of the GPPZ metric 
follows from the logarithmic terms in the coupled field equations. Combining 
with the purely gravitational result h^j, defined in ( |5.38|) , that is both 
transverse and traceless, one gets 


^l(4)y — h( 4 )ij + 1 9 -Rp’0(O) q Vj0(o) Vj0(o) T ^ (^0(0)) 9(0)ij 


12 3 

+ g0(O)VjVj0(o) + — 0(O)Do0(O)5 l (O)p- 


12 


1 


— /i( 4 )jj + -X0 + -5 , (o)p(0(o)Do0(o) + g-R0(o)) (5.58) 

where h t] is given by (|5.38|) and T.fj is the stress energy tensor derived from 
the action 

50 = f d 4 x^g^A sca i , (5.59) 
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where A sca i is the matter conformal anomaly given in (|5.66| ). Furthermore 
/r 2 ( 4 )ij = 0 in the GPPZ case. 

Using the asymptotic solution one can calculate the on-shell regularized 
action as in (|5.40| ) 


Sreg = / d x^g^) + — 0 (o) + loge ( — (Rij[9(o)]R lJ [9(o)} ~ ^R[9 (o)}) 
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-(0(O)Do0(O) + qR[9(o)}4>(o ))) + &( 


(5.60) 


To obtain the counterterm action we need to rewrite the divergences in terms 
of the induced fields 'jij and $(x, e). This amounts to formally inverting the 
asymptotic series expressing 7 ^ and <h(x, e) in terms of g( 0 ) and 0 (o) to obtain 
the sources in terms of the induced fields. Inserting this result in ( |5.60| ) (and 
including an overall minus) yields 


S t± = 




( \ - + ^ 2 (^, e) + e ) 


(5.61) 


loge 


l(/y 7 m7] - ^W 2 ) + i[$(x,e)n 7 $(x,e) + i/bh 2 (*,e)] 


The quartic term in <F, which is actually finite, is introduced to preserve 
supersymmetry [|T5] . 

We can now proceed to compute the 1-point functions. For the scalar 
operator we obtain 


(O*) 


1 ^S ren / 1 SS suh \ 

™ \ e 3 / 2 ^5$(x,e)) 


For the stress energy tensor we obtain 


-2(0(2) + 0(2)) + g0( 0 ) 


(5.62) 


(Tij) - 9(4)ij + -[Tr gf 2 ^ - (Tr g( 2 )) 2 ] 9(o)ij ~ ^9(2))^ + ^ 9 ( 2 )ij Tr g^) 

+ [0(o)(0(2) - 2 ^( 2 )) ~ -(<90 (O ))“]^(o)iy + -di4>( 0 )dj(j)( 0 ) + -(Z)“ + T^) (5.63) 


where T-‘- and 7)^ are the stress tensors of the conformally invariant trace 
anomalies. Their contribution is scheme dependent. 
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Exactly as in our discussion of the stress energy tensor in the Coulomb 
branch flow, one can use the asymptotic solution of the bulk held equations 
to compute the divergence and trace of the stress energy tensor, 


= -(O)V,<#.(„, (5.64) 

(77) = -0 ( o,(O) + - V) - h(V0 (o) ) 2 - jR^} (5.65) 


These are the expected Ward identities. The second term in ( |5.65|) , is the 
holographic gravitational conformal anomaly (5.48) [13] and the last term, 


■A SC ai — 9 [(V0(o)) 2 g-R0(o)] 


(5.66) 


is a conformal anomaly due to matter [34, 14, 35]. The integrated anomaly 


should itself be conformal invariant, and indeed (|5.66|) is equal to the La- 
grangian of a conformal scalar. The coefficients are also the ones dictated by 
non-renormalization theorems. 

In principle n-point functions involving T\j plus any set of additional 
operators containing T ki and 0$, can be computed from (|5.63|). However, 


the key information is contained in the generically non-local dependence of 
g( 4 )ij and 0( 2 ) on the sources. To first order in sources this information can 
be obtained from the linear fluctuation equations, as was done in [|15j. To 
higher order, i.e. for 3-point correlators and beyond, one needs either a 
fully non-linear solution of the held equations or a perturbative formulation 
via Witten diagrams in which linear fluctuations play the role of bulk-to- 
boundary propagators. 

Finally we note that supersymmetry requires {T tJ ) = 0 when sources are 
set to their values in the domain wall background. It is this physical re¬ 
quirement that is satisfied because of the inclusion of the Enite ‘counterterm’ 
quartic in <3>. 


6 Vector sector 


6.1 CB Flow 

In this section we study gauge fields in the fixed domain-wall background 


(|5.25|). The physics whose holographic description we seek is that of the 
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conserved currents of a spontaneously broken symmetry and their accompa¬ 
nying Goldstone bosons. The components of the operator Tr X 2 probe this 
physics. With reference to the decomposition of the 20' representation of 
SO(6) given in Sec. [2.1| , the (1,1) 0 component has non-vanishing VEV. ft is 
dual to the bulk scalar $ discussed in earlier sections on the Coulomb branch 
flow. The 8 components of TrX 2 in (2, 2) + © (2, 2)_ are interpolating fields 
for the Goldstone bosons. The bulk duals of all 20' components of TrX 2 
can be packaged as a traceless symmetric M tJ with the dynamics of a gauged 
non-linear cr-model [|l^]. The Goldstone components of M V} are related to 


the “angular” holds used in our work by (12) of ]T9|. 

We thus concentrate on the coset sector of the Coulomb branch. In our 
linearized treatment it is sufficient to use a single vector held A M (and ac¬ 
companying phase a) to describe any of the 8 vectors in this sector. It is 
the longitudinal and radial components of A p that are relevant to the broken 
symmetry physics. However it is artihcial to separate them completely from 
the transverse components and we include results on these in coset, SO(2) 
and 50(4) sectors. Final transverse correlators agree with [|T9| . 

The hrst step in the analysis is the near boundary solution of the vector 
equations of motion from the action ( |3.1| ). The functions K in the gauge 
kinetic terms all approach 1 near the boundary, i.e. K = 1 + O(p), as one 
can verify using ( |3.3| ) and ( |5.24| ). The AdS masses of the unbroken vectors 
are vanishing Mj 0 ^ = Mg 0 ^ = 0, while for the coset vectors we have the 
asymptotic expansion 


m 2 = Pp (i - -Pp) + e>(p 4 ) 


(6.67) 


An uncoupled equation for B p can be obtained by combining ( |3.13| ) and 


(|3.15|) . See |1§ for details. Its most singular solution behaves as p 1 near 


the boundary. From ( p,15|) we learn that the leading behavior of Bi is loga¬ 
rithmic. We therefore postulate the following expansions 


B P = - (B ( o) + p(B( 2 ) + log pB {2 )) + •••) (6.68) 

Bi = B( 0)i + log pB(p)i + p{B (2 )i + logpH( 2 )j) + ... (6.69) 
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Inserting these expansions in the held equations one obtains^ 

B( 2) = —-D-B(o), 0 ), diB^y = —4 (£>(2) + -8(2)) + -8(o)f? 2 

-S( 2 )i = &iB(2) ~ diB(2)i = D(-B( 2 ) — - 8 ( 2 )) — ^ 8 ( 0 ) (6.70) 

where -F( 0 )p- is the held strength of B ( 0 )j. 

We have so far presented the discussion in terms of the gauge invariant 
combination B p = A p + d p o!^\. To derive correlation functions of J* and O a 
we need to rewrite the equations in terms of corresponding sources. We will 
work in the axial gauge A p = 0. In this gauge one can integrate (| 6 . 68 | ) to 
obtain a, 


a(x,p) = log pB {0) (x)+a {0) (x)+p(B {2) (x)-B {2) (x)+log pB (2) (x))+... (6.71) 

. 8 ( 0 ) is the source for the dual operator O ai and as we shall shortly see, cq 0 ), 
the “integration constant” is related to the VEV of O a . 

From ( 6.69 ) and ( 3.71 ) we obtain 

Ai — 7l(o )i + p(/l( 2 )i + logpi ( 2 )i) + ••• (6.72) 

where the coefficients satisfy 


diA {2)i = -tB 


(0), ^(2)i - 


(6.73) 


The expansion of Ai is also valid for both longitudinal and transverse com¬ 
ponents. 

From diBi = Dcqo) and (|6.71|) we fold 


DQ(o) — — 1 diA(o)i — 4.8(2) + 0-8(0) + -FB^ 


(6.74) 


10 If one extends the system of equations to include the scalar $ one finds that there is 

a back reaction whose leading singular term is 3$ = — ^/3/2 p log 2 We have shown 

that this correction may be included in the <I> counter terms already written in (5.42). The 
back reaction then does not affect the computation of current correlators, so it is ignored 
in order to simplify the present discussion. 

11 The scalar field a is normalized so as to absorb a factor of the gauge coupling which 
is g = 2 in our present conventions. 
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The near-boundary solutions are now inserted in the action cut off at 
p — e, and one finds the following divergences 

s „g = / d , xVmtoge(~F m J^ ) - CB ( 2 0) ) (6.75) 

The counterterms are given by 

s a = j (6-76) 

We have expressed the gauge invariant quantity B( o) in terms of the cr-model 
Goldstone field m(x, e) = E A ea(x, e) at the cutoff. This is not invariant 
under the residual gauge transformation 6a = <5ct;(o), but this is expected 
as discussed below. 

We can now derive the 1-point function (in the presence of sources). By 
varying the renormalized action, with both 5Ai and 6m = e(loge48( 0 )-|-<5a;(o)) 
we obtain 

6S ien = -2 J d A x y /g^([6A i0)i + di(6a^) + loge6B w )]F pi 

—(5a(o) + loge5i?(o))f' 4 ^—b -SAfo^djF^ji logej (6.77) 

— — 2 j d A Xy^g[ o) {6A(Q)i{B(2)i + B^)i — 9*-8(2)) — ^ 4 55(o)Q((o)) 


where the first line comes from the variation of the regularized bulk action 
and the second from the variation of the counterterms. The last line was 
obtained by partial integration of the derivative of the variation 6m and use 
of ( |6.70| ). Notice that all logarithmic divergences cancel. Furthermore, the 
identification of A^i and B ( 0 ) as sources and cqo) as a VEV is now clear. 
Since the operator dual to 8 (o) has dimension 2 , the source should have 
dimension 2 too, so the dimensionally correct source is a(o) = E 2 B (op 

Since we are not interested in correlation functions of the stress energy 
tensor we can set (y'(o)p- = Sij- Going to momentum space, we rewrite ( |6.77| ) 
as 




J d A p [6A m [K ij (-2B { 2 )j - ^p 2 A Wj ) - 2i^fcp 0) \ 

+6a (0) £ 2 [-2i^A m - —(-4B (2) - ^-a (0) + —«(o))]^ (6-78) 
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where we have used ( |5.70| ) and (|6.74| ) and the dehnition of a(o). Note 
that this expression is not invariant under the residual gauge transforma¬ 
tion 5A( 0 )i = ipiSX, 5(1(0) = 0, but one obtains SS ien = —2 / d 4 x£ 2 a^). The 
failure of gauge invariance is related to the appearance of a Goldstone pole. 
A similar phenomenon occurred in the treatment of the gravitational sector 
for the Coulomb branch. (See (5.7) of |T5[. The gauge transformation there 
is the linearization of (2.17) for dg^j.). 

We can read the transverse, longitudinal and scalar 1-point functions from 

MM, 


~ 25(2 )j 2 ^ ^(°b') 

= — 2d(o) 


(6.79) 

(6.80) 


(O a ) = -2 i V -A 2 A m - =^( 




P 


P 


-45(2) - ^-O(o) + g a (o})) (6.81) 


It follows that the 2-point functions of the longitudinal currents is zero and 
the one of the current with the operator O a is trivial up to the expected 
0 -mass pole, 

(6.82) 


( Ji(p)O a (-p )) = 2 i£ 2 ^ 

pZ 


This correlator is real in .x-space, as it should. The non-trivial 2-point func¬ 
tions of the transverse current correlators and the scalar operator will be 


presented in section 8.1. 


6.2 GPPZ Flow 

The physics of this sector is that of the 5(1)^ current of an Af = 1 decom¬ 
position of J\f = 4 SYM theory. Current conservation is spoiled by the mass 
deformation, so we expect to find that 

V l Ji = q(30* (6.83) 


where O ® is a U(1)r rotation of the mass operator 0$, and q = —2/3 is the 
5-charge of O^. These operators are both in the anomaly supermultiplet of 
the mass-deformed theory, so we should obtain f3 = — \/3 as in the gravity 
sector 


15 
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The Lagrangian of this sector was derived from J\f = 8 , D = 5 gauged 
supergravity in |18| and presented both in the non-linear Stiickelberg for¬ 
malism (see (124) of |]TB[) and linearized (see eq.(42) of [|TSJ ). In our present 
notation this Lagrangian reads 


1 „ 3 » 1 

VG 4 ^ + 2 


(cLT ) 2 + 3 skill 2 


" 24> N 


- -d^af 


+ !/($) (6.84) 


where K = 3 and the vector mass is defined in (|3.2|) . Notice that we kept the 
normalization of jjHJ for the Stiickelberg field. The Stiickelberg formalism is 
convenient to solve the equations of motion, but it somewhat obscures the 
physics. To bring out the physics we first linearize near the boundary where 
<f> —> 0 , obtaining 


1 £ = -F 2 + - 
VG 4 ^2 


(^$) 2 + 4<h 2 (^--^«) 2 + V (<E>) 


(6.85) 


The change of variables described in Sec. 6.2 of [18|], which at the linearized 
level (and in our present notation) reads 

<p = $ cos(a) '0 = $sin(oi) ( 6 . 86 ) 

leads to the quadratic Lagrangian 

= \ F l» + \[( d M + 2A + i 9 ^ ~ 2A nV) 2 ] (6-87) 

The quadratic potential from (|2.1C| ) must be added. It is gauge invariant. 

A rescaling of A M is needed to bring ( |6.87| ) to a form in which the U( 1)r 
generator T rt is conventionally normalized. In |TE[ the unusual normalization 
tr 15 (Tx) = 12 in the adjoint of 577(4) was used, whereas fr 15 (T|) = 4/3 is 
conventional in J\f = 1 field theory. We therefore scale A M — * A^j 3 giving 
the Lagrangian 

+ \[{ 9 ^ + \ A ^f + ( 9 ^ - \ A M 2 } ( 6 - 88 ) 

With this form of the vector Lagrangian, the source of the standard U(1)r 
current is the new A^. The field ^ approaches the boundary at the rate 
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with = \/3 given by the background value in fl5.50| ). Thus the 
source of the operator Cf is \/3a. 

We now construct the near-boundary solution of the field equations for 
the Stuckelberg field B p = A p — f<f a. The vector equations are obtained 
from ( |6.84| ) after the rescaling. The appropriate expansions are 


Bi{x,p) = 5 (0)i + p5 (2 ), + p log p5 (2 ), + ... (6.89) 

B P (x, p) = 5 (0 ) + logp5 (0 ) + p(B( 2 ) + logpS (2 )) + ... (6.90) 


The transverse components of 5* satisfy an uncoupled equation of motion 
while longitudinal and radial components are coupled. The expansions are 
substituted and one can solve order by order in p to obtain 

B(2 ) = -^(n - 4)5(0), B(2) = (^n - !)^(0) - ^(n - 4)5(0), 

V*5( 0 )i = —45(o), V*5( 2 )i = D5(o) — 45( 2 ), (6.91) 

B(2)i = (1 - -□)5( 0 )j. 


We will work in the axial gauge A p 
(|6.90|) to obtain a(x,p), 


a{x,p) = -- 


O(o) _ B( o) + p(5(o) 


0. In this gauge one can integrate 
^5(2)) + plogp5 (0 ) + ...j (6.92) 


where cq 0 ) is an integration constant. Let us call a( 0 ) = «(o) — 5(o)- One can 
now obtain the expansion of A,(x, p) but we do not report this formula here 
since the explicit expressions for the expansion of 5,(x, p) in (|6.91| ) will be 
sufficient for our purposes. 

Substituting the near-boundary solution in the action (|6.84|) we compute 
the regulated on-shell action whose counterterms are 




1 

6 


d A x V7 l°ge 


\l ik J i l F vj F kl + iM 2 (e) 7 ^5 ? ;5 i 


As usual Srer, = Hill, 


■'ren >0[‘S'reg T “Set]- 

The variation of the renormalized action is 


(6.93) 


£Sren — —- J ffpof 5( 0 ), [5(2), + 5(2), ~ 9,5( 0 )] 


(6.94) 


i[B(2)i + 5(2), — 9,5(0)] — 65, 


(o) 


<ffi(o)) 
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where we have used 5( 0 )j = A ( 0 )i — |^a(o) and integrated by parts. Thus we 
hnally obtain 


( J %) — — g [B(2)i + B(2)i ~ 9iB(o)] (6.95) 

(O a ) = 25(0) (6.96) 


They satisfy the Ward identity 


v ‘«) = 0 ) = \(O a ) = |^3 ( 0 „) = 


(6.97) 


where in the first equality we used the (|6.95|) and (|6.91|) , in the second 
in the third the relation between the sources of O a and of Cfy derived earlier 
and in the last one the value of the /3 function. 


We have already derived the Ward identity (|6.97|) using the Stuckelberg 
gauge invariance in section 4.5, but it will be instructive to rederive it starting 
from ( |6.88|) . The Lagrangian in ( |6.88 ) is uniquely fixed by gauge invariance 
and the requirement that in the UV limit the correlators approach that of 
J\f = 4 SYM with standard normalizations. The charges of (j) and under 
U(1)r are opposite to that of the operators they couple to. The Lagrangian 
(|6.87|) is invariant under the following gauge transformations 


<V = ~3 A V5 


H = 


8 An = d fl \ 


(6.98) 


Evaluating on the background we get that the sources transform as 


(fy(o) — 0, (fyfy) — -A<E>b, 


5 = <fyA 


(6.99) 


where <p( 0 ) and ^( 0 ) are the sources for 0$ and Cfy, respectively, and <h_B is 
the background value of $. This then leads directly to the Ward identity 


v*«) = 


( 6 . 100 ) 


This derivation should be contrasted with the one of the relation (Tfy = 
/ 3(0 $) presented in ]T5[. There Weyl invariance implied (T-) = — cfyofyO#) + 


A, which upon evaluation on the background solution lead to the desired 
/3-function relation. 
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7 Linearized analysis around domain-walls 


As remarked several times, the near boundary analysis does not fix certain 
asymptotic coefficients, associated with operator vevs, in terms of the corre¬ 
sponding sources. One needs a solution of the field equations which is valid 
beyond the asymptotic region of small p. For the purpose of computing 2- 
point functions a linearized solution around a given background is enough. 
We have already used this strategy in the gravitational sector in [ 151, and 
we now extend it to vector fields. We will thus recall from [TS] the relevant 
linear equations and their solutions. 

As shown in [II 8| , for supersymmetric flows, transverse modes of the metric 


and vector fields Q can be expressed in terms of an auxiliary “massless scalar 
field” /(r). For transverse traceless tensor fluctuations one has 


x) = e i:j (p)e ip x f p (r) (7.1) 

while for transverse vectors one finds 

aj (r, x) = Vi (p) e ipx K~ l / 2 e 2A -y- (7.2) 

dr 

Longitudinal and radial modes are less universal and must be studied on a 
case by case basis. 


7.1 CB flow 


With v as radial variable, the “massless” scalar field equation has solution 


f p (v) = v a F(a, a; 2 a + 2; v ) 


(7.3) 


where 

_ + (7 - 4) 

12 Let us mention that by rescaling the vectors so that their kinetic term is canonically 
normalized, one obtains that all transverse vector fluctuations have the common mass fl8|| 


». 2 1 0 2 K 

Mitt = — 

e// 2 K 


A! 


d r K 

~K~ 


1 fd r I<\ 
4 V“fTj 


M 2 

~K 


= -28 2 A. 


We have checked that this relation is satisfied by the graviphoton in the GPPZ flow and 
by any vector field in all known supersymmetric CB flows with one active scalar. 
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The transverse vector fluctuations thus read 

( v ) = v x+a ~ 1 F(a, a + 1; 2a + 2; v) (7.5) 

where A = 1 for 50 (4), A = 2 for 50(2) and A = 3/2 for the coset. The 
‘unbroken’ 50(4) x 50(2) vectors have no physical longitudinal fluctuations. 
The longitudinal components of the coset vectors can be expressed in terms 
of the radial component B r . In turn, B r can be rewritten as 

B r {r,x) = e~ AA M~ 2 C p {r)e ip ' x (7.6) 

When expressed in terms of the v variable, the scalar function C p {y) satisfies 

v'\l - v)c; + v(l - v)C r - + (1 - v)]C r = 0 (7.7) 

whose solution is the hypergeometric function 

C p {y) = v a+1 ^ 2 F(a, a + 1, 2a + 2; v) (7.8) 

with a defined in (04). Up to (irrelevant) overall constant, 

C p [y) = v(l - v)B v (7.9) 

The coordinate transformation (|5.24j) may be used to obtain B p . 


7.2 GPPZ flow 

When expressed in terms of the radial variable u, the “massless” scalar field 
equation admits the solution 


/ p (n) = (l-n) 2 F(2 + z|,2-*|;2;n). 

For the transverse components of the U( 1)r graviphoton one has 


b^\u) = u( 1 - u)F (2 + 2- 3; u) 

after dropping an irrelevant p-dependent factor. 

For the longitudinal and radial components one finds 

'3131 


■P 


C P {u) = SB p = u f(^- + -q, - - -q- 3; u) 


where 


q = Vi-p 


(7.10) 


(7.11) 


(7.12) 

(7.13) 
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8 Correlation functions of vector currents 


In this section we will obtain results on the 2-point correlation function of 
various vector currents and their scalar partners in the CB and GPPZ flows. 
We will study this system in the fixed background approximation in which 
the back reaction on the background metric and active scalar is neglected. It 
was argued earlier that the fixed background approximation is sufficient to 
compute all correlators in this sector except those involving T tJ . 

It is actually easy to see that the 2-point function (T l3 (p) J k ( —p)) = 0, 
so that nothing is missed in the fixed background treatment. On the grav¬ 
ity side this correlator could only come from bilinear terms h^A *. in S Ten , 
and such terms are clearly absent. For the held theory of Coulomb branch 
how, (Tij(p)Jk(—p)) must be conserved in all 3 indices, and there is no ten¬ 
sor structure with this property. For the GPPZ how the tensor structure 
(Tij(p)Jk(—p)) = 7TijPkA(p 2 ) is allowed. However we trace on ij and contract 
with pk and use the operator equations of the theory we fold that A(p 2 ) is 
proportional to the 2-point function of a scalar and pseudoscalar operator 
which vanishes by parity conservation. 


8.1 CB flow 


For the CB how the physics of the ‘unbroken’ and ‘coset’ currents, is dif¬ 
ferent. The distinction between sectors is almost irrelevant for transverse 
components. Since only coset vectors have longitudinal components, we can 
treat all sectors together with a small elaboration of notation. 

The hrst step is to perform the asymptotic expansion of the hypergeo¬ 
metric functions that appear in the solutions of the linearized held equations. 
From (73) one hnds that the transverse vector fluctuation behaves as 


A — A(p)i ^1 - pf[ A - 1 + - |^(^( a + ^ ~~ + 2 + • • J 

( 8 . 1 ) 

where A = 3/2 for the coset vectors, while for the ‘unbroken’ vectors A = 0 
for SO(2) and A = 1 for 50(4). For the radial component in the coset sector, 


one can use (|7i8|) to obtain 


5(2) — O(o) 


1 

6 


^ + ^(a+l)-W 


( 8 . 2 ) 
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This information is sufficient to fix B^)i and B( o) which were the two unknown 
coefficients in ( 6.78 ). This variational expression may now be “integrated” 
to obtain the bilinear fluctuation action 


Q 


N 2 
202 


J d A p A {0)i (p)TT i:j A {0)i (-p ) \ -i 2 { A - 1 ) + y [ 0 (a + 1 ) - - 0 ( 1 )] 


p 


2£ 2 


8 x 2it o ( o) (p)—A( 0 )i(—p) + [-- + 2(0(a + 1) - 0(l))Ho)(p)a(o)(-p) (8.3) 


P 


p 


where we have reinstated the overall factor of N 2 /2ty 2 , and where S x = 0 
for the ‘unbroken’ /SO(4) x SO( 2 ) vectors reflecting the absence of physical 
longitudinal modes, while 8 x = 1 for the coset vectors. 

The negative of the second variation of S Ten wrt to sources gives the the 
desired 2-point functions. For the /SO(4) x SO(2) longitudinal modes one 
obviously gets 

(J‘(p) ^(-p ))(0 = 0 (8.4) 

since there is no mixing with scalars. For the 2-point function of current and 
the scalar operator, we can repeat ( |6.82| ), slightly generalized to read 


/V 2 nr \. 

{Mp)OJ-p)) = ^S x -2U 


(8.5) 


The 2-point function of the operators dual to the coset scalars is obtained 
from ( B-3p . It reads 


{0 a {p)0 Q (-p)) = ^ - 4(0(a + 1) - 0(1)) j ( 8 . 6 ) 

This has the expected massless pole with a residue proportional to £ 2 and 
also the correct UV behavior. 

The transverse current correlator is 

N 2 

( Ji(p)Jj{q))(t ) = ij (- 200 A - 1 ) + P 2 [0(a + 1 ) - 0 ( 1 )]) . (8.7) 

It has the correct scaling in the UV but it also has a massless pole except 
for A = 1, i.e. for the 50(4) currents. The result agrees with |19|. Although 
the supergravity solution (’continuous’ disk distribution of D3-branes) seems 
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to preserve an 50(4) x 50(2) subgroup of 50(6), no discrete disk distri¬ 
bution can preserve more than 50(4). The field-theory counterpart of this 
phenomenon is the absence of any locus in the Coulomb branch of J\f = 4 
SYM with any finite number of colors N on which 50(4) x 50(2) is pre¬ 
served. This is discussed in the Appendix. We thus interpret the massless 
poles as Goldstone poles in both the coset and 50(2) sectors. The latter 
indicates that supergravity implicitly agrees with held theory, in accord with 

10 . 


the discussion in the Addendum to 


8.2 GPPZ flow 


In this section we combine results of near boundary analysis from section |T2 


with the fluctuations given in section o to obtain correlation functions of 
the U(1)r current and the scalar operator O a . 

It is again useful to ‘integrate’ the first order variation 5S ien in ( J6.94|) to 
obtain a quadratic action from which the 2-point correlators may be immedi¬ 
ately read. We first determine the non-local relation between 5(2)*+ 5(2)* and 
A( 0 )i from the asymptotics of the hypergeometric functions given in section 
|as«-> 1. For the transverse components (|7.11|) gives 

Bi = A (0 )j(l - p[(l + )[(A - -) - logp)] + 1] + •••) (8.8) 

where K = -0(3) + -0(1) — 0(2 + ip/2) — 0(2 — ip/2). One then finds 

2 2 1 
5(2)* = Ao)i(l + ^-) 5( 2 )i = —A(0)*[(1 + ^-)(A - -) + 1] (8.9) 

Similarly, for the radial component, (|7.12| ) gives 

5(o) = — 5(o) J, (8.10) 

where J = 20(1) — 0(3/2 + q/2) —0(3/2 — q/2) and from the near-boundary 
analysis we know that 

V*5( 0 )j = —45( 0 ) (8.11) 

The renormalized action to quadratic order in the sources thus reads (in 
momentum space) 


Q 


N 2 
12?r 2 


d a P Uo)XW)i[-(l + ~\) + J 
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1 


( 8 . 12 ) 


{PiA\ o) + ~ip 2 a {0) ) 2 J 


2 p 2 ^ (0) ' 3 


The fact that the the longitudinal part of the and the a(o) appear as 
a total square is a consequence of the beta function operator relation ( |8.17|) 
between the divergence of J* and O a . We have seen an analogous phenomenon 
in the graviton-scalar section in [|15[| . 

The transverse 2-point function reads 

(■((1 + j)(K ~ 1) - £) (8.13) 


It has poles at p 2 = —4 (n + 2 ) 2 (we are using Euclidean signature) with 
n = 0, 1 ,as expected, but also a disturbing massless pole whose residue 
is —(N 2 /6tt 2 )(K(p=0) — 1/2) = ( N 2 /6tt 2 ). Happily, as we now show, the 
longitudinal 2 -point function also contains a massless pole and the two con¬ 
tributions cancel each other! 

The remaining correlators are 

(Mp)JA-p)) m = ( 8 - 14 ) 

( Ji(p)O a (-p )) = (8.15) 

N 2 „ _ 

{O a (p)O a (-p)) = J ( 8 - 16 ) 

The residue of the zero mass pole in (Ji(p)Jj(—p)){i) is N 2 /(6tt 2 ), and indeed 
the zero mass poles cancel. 

These correlation functions are consistent with the operator relation 

V* J* = \0 a = (8.17) 

where (3 = —y/3 is the same beta function found in the graviton-scalar sector 
in |H|. The current Ward identity ( p.97| ) is related by supersymmetry to the 
trace Ward identity, T/ = (30 $. 


Acknowledgments 


We would like to thank M. Berg for a discussion of the results presented in 


21 ] prior to publication and for finding some typographical errors in an early 


41 













draft of this paper. Part of the present work was done during the summer 
workshop “Physics in the Pyrenees: Strings, Branes and Field Theory 1 ' held 
in Benasque, Spain. We thank the local organizers and other participants 
for the stimulating atmosphere there. The collaboration of M.B. and D.Z.F. 
was facilitated by the INFN-MIT ‘Bruno Rossi’ exchange program. Much 
of this work was carried out while M.B. was visiting D.A.M.T.P. at Cam¬ 
bridge University. M.B. would like to thank the colleagues at D.A.M.T.P. 
and in particular Michael Green for their kind hospitality and for financial 
support through a PPARC grant. The work of M.B. was supported in part by 
the EEC contracts ffPRN-CT-2000-00122 and ffPRN-CT-2000-00148, and 
in part by the INTAS contract 99-1-590. DZF’s research is supported by US 
NSF grant PHY-0096515. KS is supported in part by the National Science 
Foundation grant PHY-9802484. 

A Appendix 

In this Appendix we use simple group-theoretic methods to prove a lemma 
concerning the Coulomb branch of J\f = 4 SYM theory with any finite di¬ 
mensional gauge group G. 

Lemma: The Coulomb branch does not contain configurations where the 
50(6) flavor symmetry is broken to a non-trivial direct product SO{n ) x 
50(6 — n ) with 2 < n < 4. 

The proof assumes the validity of a classical description of the moduli 
space in terms of the VEV of the elementary scalars X 1 = X r*T a . The VEV 
(X*) is in the Cartan subalgebra (A"*) = 

Without loss of generality we can consider the subgroup 50(2) C 50(n) 
which rotates in the i — 1,2 plane. This acts as a complex rotation of the 
sum Z = X 1 + iX 2 , and preserving 50(n) symmetry requires 

e iip Z — gZg- 1 (A.18) 

where g is a gauge transformation depending on ip. For infinitesimal 8(p this 
requires a Lie algebra element A such that 

6(pZ=[A,Z\ (A. 19) 

We can write A = zb + ) + A^) + A^ H \ where zb + ) belongs to the positive 
root space, A^> is its adjoint in the negative root space and is in the 
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Cartan subalgebra. Thus 


8ipZ = [A, Z] = S (+) - B { - ] (A.20) 

where RU) are in the positive (negative) root spaces and are linear in Z. 
However the last equation is impossible to satisfy unless Z = 0. 

Thus for any SO (n)-invariant vacuum, the non-vanishing components of 
X 1 are restricted to % — n + 1, .. 6 . the same argument may now be applied 
to these components which also must vanish. 

We would now like to rephrase the above discussion in terms of gauge 
invariant composite operators. In this approach the question is whether it 
is possible to pick VEV’s of the elementary fields X 1 such that the only 
composite operators that acquire a VEV are singlets under SO(n) x SO (6 — 
n ). The composite operators that admit a supergravity description are the so- 
called chiral primary operators (CPO). They are of the form Tr ( X k ), i.e. k- 
fold symmetric traceless products of the fundamental of AO( 6 ), that belong 
to the representation with Dynkin labels [0, k,0]. Under the decomposition 
of SO( 6 ) R-symmetry into SO(n) x SO (6 — n) the operators with k even 
contain singlets under SO(n) x SO (6 — n) but the ones with k odd contain 
singlets of each factor only but not of the product. So in order to preserve 
SO(n) x SO (6 — n) the VEV’s of Vs should be such that (i) for even k only 
the singlet operator gets a VEV, (ii) the VEV of all CPOs with odd k are 
zero. We will now show that it is not possible to pick VEV’s such that these 
conditions are satisfied. 

Let us consider the case of a U(N) gauge group. We will comment on the 
generalization to SU(N) afterwards. The lowest CPO’s in the 20' of SO( 6 ) 
contain one singlet of SO(n) x SO (6 — n), say Tr ((6 — n)[Xf + ... A -2 ] — 
n[A " 2 +1 + ... X|]). In order for the other non-singlet components to vanish 
one has to choose the scalar VEV’s so that the six iV-component vectors 
Xi be orthogonal to one another and normalized so that |Ab | 2 = ... = 
|X n | 2 7 ^ |AV+i | 2 = ... = AT 6 1 2 . Let us now turn to the cubic Casimir 
Tr (XiXjXk — ■ ■ •) in the 50 of SO( 6 ). It does not contain singlets of the 
product subgroup SO(n) x SO (6 — n), so all VEV’s of the CPO in the 50 
should vanish. The only way to achieve this is to take X t along the generator 
Hi and Xi = ... — X n =/= X n+1 = ... = X e . The next CPO in the 105 
of SO( 6 ) is associated to the quartic Casimir Tr (X^X^X^Xi — ...). Let 
us consider the complex field Z = X 1 + LA 6 . It is easy to check that the 
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operator Tr [ Z 4 \ is not a singlet of the product subgroup, and one cannot 
further restrict the VEVs of X 1 (without setting all of them to zero) such 
that it vanishes. We thus conclude that R-symmetry cannot be broken to 
SO(n) x 50(6 —n), and at least one of the two factors, say 50(6 —n), should 
also be broken. In the case of SU(N) the only additional complication is that 
the tracelessness condition Tr Hi = 0 should be taken into account. This 
imposes a further constraint among the VEV’s that is however compatible 
with the previous choices. 
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